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The thermodynamic and dynamical properties of an Ising model with both short range and long 
range, mean field like, interactions are studied within the microcanonical ensemble. It is found that 
the relaxation time of thermodynamically unstable states diverges logarithmically with system size. 
This is in contrast with the case of short range interactions where this time is finite. Moreover, at 
sufficiently low energies, gaps in the magnetization interval may develop to which no microscopic 
configuration corresponds. As a result, in local microcanonical dynamics the system cannot move 
across the gap, leading to breaking of ergodicity even in finite systems. These are general features 
of systems with long range interactions and are expected to be valid even when the interaction is 
slowly decaying with distance. 
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Systems with long range interactions are rather com- 
mon in nature jlj . In such systems the inter-particle po- 
tential decays at large distance r as 1/r" with a < d, 
in dimension d. Examples include magnets with dipolar 
interactions, wave-particle interactions gravitational 
forces 13], and Coulomb forces in globally charged sys- 
tems 4J. It is well known that in such systems the vari- 
ous statistical mechanical ensembles need not be equiva- 
lent . For example whereas the canonical ensemble spe- 
cific heat is always non-negative, it may become negative 
within the microcanonical ensemble when long range in- 
teractions are present Q . It has recently been suggested 
that whenever the canonical ensemble exhibits a first or- 
der phase transition the canonical and the microcanoni- 
cal phase diagrams may be different 0. This has been 
demonstrated by a detailed study of a spin-1 Ising model 
with long range, mean field like, interactions. While the 
thermodynamic behavior of such models is fairly well 
understood, their dynamics, and the approach to equi- 
librium, has not been investigated in detail so far [8|. 
The aim of the present Letter is to identify some general 
dynamical features which distinguish systems with long 
range interactions from those with short range ones. 

One characteristic of systems with short range inter- 
actions is that the domain in the space of extensive ther- 
modynamic variables X over which the model is defined 
is convex for sufficiently large number of particles. Here 
the components of the vector X are variables like the 
energy, volume and possibly other extensive parameters 
corresponding to the system under study. Convexity is a 
direct result of additivity. By combining two appropri- 
ately weighted sub-systems with extensive variables Xi 
and X2 , any intermediate value of X between Xi and X2 
may be realized. On the other hand systems with long 
range interactions are not additive, and thus intermedi- 
ate values of the extensive variables are not necessarily 



accessible. This feature has a profound consequence on 
the dynamics of systems with long range interactions. 
Gaps may open up in the space of extensive variables and 
lead to breaking of ergodicity under local microcanoni- 
cal dynamics of such systems. An example of such gaps 
in a class of anisotropic XY models has recently been 
discussed in Q. 

Another interesting feature of systems with long range 
interactions is their relaxation time. It is well known 
that the relaxation time of metastable states grows ex- 
ponentially with the system size On the other hand 
the relaxation processes of thermodynamically unstable 
states are not well understood. In systems with short 
range interactions the relaxation takes place on a finite 
time scale. However previous studies of a mean field XY 
model suggest that this relaxation time diverges with the 
system size 0- 

In the present Letter we study some of the issues dis- 
cussed above by considering a spin-1/2 Ising model with 
both long range, mean field like, and short range near- 
est neighbor interactions on a ring geometry 0, [3| • We 
study its thermodynamic and dynamical behavior in both 
the canonical and microcanonical ensembles. It is found 
that the two ensembles result in different phase diagrams 
as was observed in other models 7]. Moreover, we find 
that for sufficiently low energy, gaps open up in the mag- 
netization interval — 1 < m < 1, to which no microscopic 
configuration corresponds. Thus the phase space breaks 
into disconnected regions. Within a local microcanonical 
dynamics the system is trapped in one of these regions, 
leading to a breakdown of ergodicity even in finite sys- 
tems. In studying the relaxation time of thermodynam- 
ically unstable states, corresponding to local minima of 
the entropy, we find that unlike the case of short range 
interactions where this time is finite, here it diverges log- 
arithmically with the system size. We provide a simple 
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FIG. 1: The canonical and microcanonical {K,T) phase di- 
agram. In the canonical ensemble the large K transition is 
continuous (bold solid line) down to the trictitical point CTP 
where it becomes first order (dashed line). In the micro- 
canonical ensemble the continuous transition coincides with 
the canonical one at large K (bold line). It persists at lower 
K (dotted line) down to the tricritical point MTP where it 
turns first order, with a branching of the transition line (solid 
lines). The region between these two lines (shaded area) is 
not accessible. 



explanation for this behavior by analyzing the dynamics 
of the system in terms of a Langevin equation. 

We start by considering an Ising model defined on a 
ring of N spins Si = ±1 with long and short range inter- 
actions. The Hamiltonian is given by 

1=1 \j=i / 

where J > is a ferromagnetic, long range, mean field 
like coupling, and K is a nearest neighbor coupling which 
may be of either sign. The canonical phase diagram of 
this model has been derived in the past 0, . It has 
been observed that in the {K, T) plane, where T is the 
temperature, the model exhibits a phase transition line 
separating a disordered phase from a ferromagnetic one 
(see Fig. The transition is continuous for large K, 
where it is given by (3 = e~^^ . Here /3 = l/T and 
J = 1 is assumed for simplicity. Throughout this work we 
take /cb = 1 for the Boltzmann constant. The transition 
becomes first order below a tricritical point located at 
an antiferromagnetic coupling Kqtp = — ln3/2\/3 ~ 
-0.317. 

We now analyze the model within the microcanonical 
ensemble. Let U — —i^ {SiSt+i — 1) be the number 
of antiferromagnetic bonds in a given configuration char- 
acterized by iV+ up spins and N- down spins. Simple 
counting yields that the number of microscopic configu- 
rations corresponding to {N^, 7V_, U) is given, to leading 
order in A^, by 

f^(A^.,iV_,C/)=(^+)(^/2) • 

Expressing 7V_|_ and iV_ in terms of N and the magnetiza- 
tion M = N+ - N^, and denoting m = M/N, u = U/N 



and the energy per spin e = E/N, one finds that the 
entropy per spin, s(e, m) = In 17 is given in the ther- 
modynamic limit by 

s(e,m) — — (1 + m) ln(l + m) + — (1 — m) ln(l — m) 

— wlnu — — (1 + m — u) ln(l + m — u) 

— -(1 — m — u) ln(l — m — u) , (3) 
where u satisfies 

-^m'^ + Ku. (4) 

By maximizing s(e,m) with respect to m one obtains 
both the spontaneous magnetization ms(e) and the en- 
tropy s(e) = s(e, ms(e)) of the system for a given energy. 
The temperature, and thus the caloric curve, is given by 
l/T = ds{e)/de. A straightforward analysis of ||2Jl yields 
the microcanonical {K, T) phase diagram of the model 
(Fig. where it is also compared with the canonical 
one. It is found that the model exhibits a critical line 
given by the same expression as that of the canonical 
ensemble. However this line extends beyond the canoni- 
cal tricritical point, reaching a microcanonical tricritical 
point at Kmtp — —0.359, which is computed analyti- 
cally. For K < Kmtp the transition becomes first order 
and it is characterized by a discontinuity in the temper- 
ature. The transition is thus represented by two lines 
in the {K, T) plane corresponding to the two tempera- 
tures at the transition point. These lines are obtained by 
numerically maximizing the entropy 

We now consider the dynamics of the model. This is 
done by using the microcanonical Monte Carlo dynamics 
suggested by Creutz In this dynamics one samples 
the microstates of the system with energy E — Ed with 
Ed > by attempting random single spin flips. One 
can show that, to leading order in the system size, the 
distribution of the energy Ed takes the form 

P{Ed) - e~^°/^ . (5) 

Thus, measuring this distribution, yields the temperature 
T which corresponds to the energy E. 

In applying this dynamics to our model one should 
note that to next order in N the energy distribution is 
given by P[Ed) e^p {-Ed/T - El/2CvT^), where 
Cv = 0{N) is the specific heat. In extensive systems 
with short range interactions the specific heat is posi- 
tive and thus the correction term does not modify the 
distribution function for large N. However in our sys- 
tem Cv can be negative in a certain range of K and E. 
It even becomes arbitrarily small near the MTP point, 
making the next to leading term in the expansion arbi- 
trarily large. This could in principle qualitatively modify 
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the energy distribution. However as long as the entropy 
is an increasing function of E, namely for positive tem- 
peratures, the distribution function ((SJ) is valid in the 
thermodynamic limit. This point is verified by our nu- 
merical studies, where an exponential distribution of the 
energy Ed is clearly observed. 

We now address the issue of the accessible magneti- 
zation intervals in this model. We find that in certain 
regions in the {K, E) plane, the magnetization m can- 
not assume any value in the interval 1). There exist 
gaps in this interval to which no microscopic configura- 
tion could be associated. To see this, we take for simplic- 
ity the case > . It is evident that the local energy 
U satisfies < C/ < 2N^. The upper bound is achieved 
in microscopic configurations where all down spins are 
isolated. This implies that < u < 1 — to. Combining 
this with one finds that for positive to the accessible 
states have to satisfy 

TO < V— 2e , m > m+ , to < to^ 

with m± =-K± ^K^ - 2(e - K) . (6) 

Similar restrictions exist for negative to. These restric- 
tions yield the accessible magnetization domain shown in 
Fig. 121 for K = —0.4. It is clear that this domain is not 
convex. Entropy curves s(m) for some typical energies 
are given in Fig. demonstrating that the number of 
accessible magnetization intervals changes from one to 
three, and then to two as the energy is lowered. 

This feature of disconnected accessible magnetization 
intervals, which is typical to systems with long range in- 
teractions, has profound implications on the dynamics. 
In particular, starting from an initial condition which lies 
within one of these intervals, local dynamics, such as the 
one applied in this work, is unable to move the system to 
a different accessible interval. Thus ergodicity is broken 
in the microcanonical dynamics even at finite N . 

To demonstrate this point we display in Fig. 0] the 
time evolution of the magnetization for two cases: one in 
which there is a single accessible magnetization interval, 
where one sees that the magnetization switches between 
the metastable m — Q state and the two stable states 
TO = ±TOs. In the other case the metastable to = state 
is disconnected from the stable ones, making the system 
unable to switch from one state to the other. Note that 
this feature is characteristic of the microcanonical dy- 
namics. Using local canonical dynamics, say Metropolis 
algorithm , would allow the system to cross the for- 
bidden region (by moving to higher energy states where 
the forbidden region diminishes), and ergodicity is re- 
stored in finite systems. However, in the thermodynamic 
limit, ergodicity would be broken even in the canonical 
ensemble, as the switching rate between the accessible 
regions decreases exponentially with N . 

We conclude this study by considering the life time 
t{N) of the TO = state, when it is not the equilib- 
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FIG. 2: Accessible region in the (m, e) plane (shaded area) 
for K — —0.4. For energies in a certain range, gaps in the 
accessible magnetization values are present. 
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FIG. 3: The s(m) curves for K = —0.4, and for typical 
energy values, demonstrating that gaps in the accessible states 
develop as the energy is lowered. 



rium state of the system. In the case where to = is a 
metastable state, corresponding to a local maximum of 
the entropy (such as in Fig. Qt) we find that the life time 
satisfies r ~ e^^** where As is the difference in entropy 
of the m — state and that of the unstable magnetic 
state corresponding to the local minimum of the entropy 
(see Fig. |5ji)- Such exponential dependence on N has 
been found in the past in Metropolis type d yna mics of 
the Ising model with mean field interactions LCl. Simi- 
lar behavior has been found in the XY model |l5| and in 
gravitational systems |l6l | when microcanonical dynamics 
was applied. 

We now turn to the case where the to = state is ther- 
modynamically unstable, where it corresponds to a local 
minimum of the entropy. In systems with short range in- 
teractions, the relaxation time of this state is finite. Here 
we find unexpectedly that the life time diverges weakly 
with N, t{N) - \og{N) (see Fig. This is to be 

compared with studies of the life time of the zero magne- 
tization state in the XY model under similar conditions 
which show that t{N) ^ with a ~ 1.7 

In order to understand this behavior we consider the 
Langevin equation corresponding to the dynamics of the 
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FIG. 4: Time evolution of the magnetization for K = —0.4 (a) 
in the ergodic region (e = —0.318) and (b) in the non-ergodic 
region (e = —0.325). The corresponding entropy curves are 
shown in the inset. 




FIG. 5: Relaxation time of the m = state when this state is 
(a) a local maximum [K — —0.4, e — —0.318) and (b) a local 
minimum of the entropy {K — —0.25, e — —0.2). 



system. The evolution of m is given by 

dm ds 
dt dm 



where ^(t) is the usual white noise term. The diffusion 
constant D scales as ~ 1 /N. This can be easily seen by 
considering the non-interacting case in which the magne- 
tization evolves by pure diffusion where the diffusion con- 
stant is known to scale in this form. Taking s(m) ~ am^ 
with a > 0, making the m — state thermodynamically 
unstable, and neglecting higher order terms, the distri- 
bution function of the magnetization, P(m,t), may be 
calculated. This is done by solving the Fokker-Planck 
equation corresponding to (j?)). With the initial condi- 
tion for the probability distribution P{m,0) — S {ml, the 
large time asymptotic distribution is found to be 

P(m,t) 



m 



exp 



D 



(8) 



This is a Gaussian distribution whose width grows with 
time. The relaxation time corresponds to the width 
reaching a value of 0(1), yielding r ^ — ln_D ~ InA^. 
Similar analysis and simulations can be carried out for 
the canonical ensemble yielding the same divergence. 



In summary, some general features of the dynamical 
behavior of systems with long range interactions were 
studied using the microcanonical local dynamics of an 
Ising model. Properties like gaps in the accessible mag- 
netization interval and breaking of ergodicity in finite 
systems have been demonstrated. We also find that the 
relaxation time of an unstable to = state, correspond- 
ing to a local minimum of the entropy, is not finite but 
rather diverges logarithmically with N. We expect these 
phenomena to appear in other systems with long range 
interactions which are not necessarily of mean field type. 
This study is thus of relevance to a wide class of physi- 
cal systems, such as dipolar systems, self gravitating and 
Coulomb systems, and interacting wave-particle systems. 
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